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ABSTRACT 


In a recent review article on the semi-empirical 
all-valence-electron method, Klopman wrote of ".... the 
painful dilemma of having to choose between a method 
that gives good heats of formation and poor bond 
distances, or poor heats of formation and good bond 
distances. ... Lts solution, if it can be found, 
would probably be the most important contribution in 
the field of semi-empirical calculation of en 
molecules." 

We have set out to meet the challenge. In the 
present thesis, the first stage of our effort will be 
described and moderately successful numerical results 
will be reported. We start with a critical review of 
the past developments of the subject and the fundamental 
theoretical background and then proceed on to the 
description of a new approach to the subject. 

A careful analysis of the total energy of molecular 
systems has led us to the point of view that a proper 
balancing mechanism between contributions from the one- 
electron and two-electron parts of the total energy is 
Vital lysimpoLtanteln Order co Sntead good heats of 
formation and good bond distances. Keeping this goal in 


mind, we have introduced a method. at the CNDO-level which 
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has some new features. A new formula for estimating 
coulomb integrals has been proposed here, This formula 
introduces two new parameters which have sufficient 
flexibility to continuously cover the rigorous formulas 
of Roothaan or the approximate formulas of Ohno and 
Mataga-Nishimoto. Another new feature is the way in 
which the one-electron part of the Hartree-Fock operator 
is parametrized. The potential of these new parameters 
for achieving the flexibility needed to accomplish the 
required delicate balance is greater than that in 
previous CNDO methods. 

rOur computer program is’ constructed at the level 
of so-called CNDO (Complete Neglect of Differential 
Overlap) level, but it is much more flexible in many 
respects than the widely distributed CNDO/2 computer 


program, 
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Chapter I 


INTRODUCTION 


Although for nearly fifty years the principles 
of quantum mechanics, which are necessary for the 
understanding of the electronic structure of molecules, 
have been known, still it has been only in recent years 
Ehategrealeadvances in abe im2tlo Calculations. OL 


(1) 


molecular electronic structure have been achieved. 
Pomc emO es LiecealeCelbecacilevementa, aD LNLuLLO 
calculations are of limited use since they are very 
expensive and require large computation time. Even if we 
spend hundreds of hours on a large computer, we can 
perhaps perform an energy-calculation of a molecule 
not very much larger than benzene and the accuracy is 
not necessarily satisfactory. The ground state 
dissociation energies of many diatomic molecules 
calculated by the Hartree-Fock method have been shown 
to be one quarter to one half of their experimental 
values. For the fluorine molecule the Hartree-Fock 
calculations do not predict any binding energy 

(-1.37 ev) whereas the experimental binding energy is 


(2) (2a) 


1.68 ev. Wahl pointed out that there might be two 
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for many molecular systems the molecular orbital wave 
function does not exhibit the proper dissociation behavior 
(i.e. it goes into neutral and ionic atomic states). 
Second, when bonding occurs the correlation energy of 

the molecule increases, and this increase tends to raise 
the electronicaceneraqy of thesmolecule. Using an optimized 


(2b) have obtained 


configuration function, Das and Wahl 
some binding energy (0.5 ev) for the fluorine molecule. 

As ab initio calculations are not possible for very 
large molecules organic and inorganic chemists use semi- 
empirical quantum chemical methods to interpret and cor- 
relate experimental results. The basic merit of the 
Semi enp rical japproachslies in its practical applicabaLiry 
to large molecular systems. Generally speaking the semi- 
empirical methods have been reasonably successful. A 
review of their development will be presented in order 
LO=Seu, Lie stege for a new formulation of san all—-valence— 
electron Semi-empirical method to be described in ‘the 
present work. 

Hund and Mulliken introduced the molecular orbital 
(MeOn) theory asvcarly as, 1928.) sin 192) Hiickei ‘3? 
applied the theory in its simple form to conjugated 
hydrocarbons. Later the theory was developed by Lennard- 


(4) (5) 


Jones, Coulson and Longuet-Higgins and others, and 


it has been successfully used to study the properties of 
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conjugated molecules, olefins and aromatic compounds. 

In Huckel naelecteuon theory wtne first, basic 
assumption is that interactions between t-electron and 
other electrons are comparatively small. w1-electrons 
are assumed to be well separated from other electrons. 
The second assumption is that the one-electron 


Hamiltonian is represented by 


TG 2 ykersteiets es h 


i=1 


efr (1) (>t) 


which means that the electrons are supposed to move 
independently in 4 potential field in which electron 
interactions are accounted for in an averaged form. Here 
the core is formed by all the nuclei and all the electrons 
other than the t-electrons. Each n-electron is treated 
as though it were in the averaged potential field due to 
the core and other 1~electrons. 

According to the assumed Hamiltonian (I-1) 
the many-body Schrodinger equation can be separated into 
a set of one-electron equations and the total wave function 


can be expressed as a simple product of MO's ss 


ery Wee I Wee) (I-2) 
Each molecular orbital is an eigenfunction of nore: 
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The total mt-electron energy is given by 


E | = 2 mies (I-4) 


where n is the occupation number of orbital Ys and can 
Deregualearer 0, 1.0 2. 

the thard basic assumption is that the MO's are 
expressed as linear combinations of atomic orbitals, i.e. 


¥, = de ee Gee) 


where the ne are the 2p. aLONLCwOLDi tals.) Therorbital 
energies and the coefficients Ce can be obtained by 


solving the secular equation 


Dei Ae Esowny te she. (I-6) 


where the matrix components are 


he, = Sida It Oa lites ae (I-7) 

Ey = 2 Wee eles? - Paro ae 
and 

aoe = 5X5) ieee (I-9) 


The secular equations are simplified by means of 


the following approximations: 





. od. WevEp-ci. Yous) tevotesiees z 


(b-T) S ih = oe ; 


(@-T) f y 2 . a 


i! 


inuMmo. x2 Se 6 


(1) Alt the comlomb imtegrals, Lana are assumed to be 


equal, i.e. Ce a en = 


(2) The resonance integrals sie are eee to be non- 
zero if atom u is bonded to atom \y, represented by 

er Ve a is set equal to zero when atom u is not 
bonded to atom vy. This approximation is known as the 


neglect of non-neighbour interactions. 


(3) All the resonance integrals for bonded carbon atoms 
have the same magnitude and are represented by the symbol 


pee Teen, Dic = IS) a 0 aye u ae Ve 


(4) All the overlap integrals 24 for u # v are neglected. 
For normalized atomic orbitals Fem ee eeyeeienesS assumed 
thatass = 6 where 6 is the Kronecker delta. MThis 

uv uv uv 


approximation is known as the neglect of overlap integrals 


or the zero overlap approximation. 


(5a) 


Coulson and Longuet-Higgins have shown that 
the total t-energy is related to n-electron charge 


density and bond order and may he expressed by 


hie 2 qs Deo Me (a0) 
u u uv uv 


where Ts and Pe are respectively charge density and bond 


order defined by 


2 
= 2) om a=) 
Fy, z «+H 
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The Huckel procedure is very Simple and ground state 
properties for conjugated molecules, especially alternant 
hydrocarbons, can be predicted well. But it does not 
give good results for ionization potentials, electron 
affinities and excitation energies. The defect in this 
simple theory is that it does not give a precise definition 
of the one-electron Hamiltonian. Electron interaction 
Homo bso motexpliciuclyetaken into account. AS a result 
there is no energy difference between a singlet and a 
triplet state. 

The next attempt to develop the n-electron 


(6) 


theory was made by Goeppert-Mayer and Sklar. Electron 


interaction was taken into account and the LCAO approxima- 


tion |equation (I-5)] was employed. The wave functions were 


made antisymmetrical with respect to exchange of electrons 
DYAMuULtap Lying. Lhemiby appropriate span functions. The 


Hamiltonian is represented by the equation 


(Gibney 
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where = is the repulsion between 1 electrons u and v 
Uv 


and ae ey ae, is the kinetic energy operator for electron uy 
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and its potential energy operator in the field of the core. 
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hice at See) + Mdealees (u) (I-14) 


OF 


Here) gh) 2 Un) + De ou) +200 * Cu) 
q#p He 


(I-15) 
where T is the kinetic energy Ihe te is the potential 


energy operator due to the ok carbon center in the core 
and Uh is the potential eer operator due to the see 
hydrogen atom in the core; * signifies a potential due to 
ey neutral atom. The atomic orbitals oe are now taken to 


be atomic eigenfunctions of the operator (T + uy 


such that 


OT ates ets el Xp) = E aycay (I-16) 


where ue an appropriate valence state ionization 
potential. The potential energy Sig tetereneienets Nats may be 


expressed as 


Cpl = ee ~fig fa) 


According to the Goeppert-Mayer and Sklar formulation 





dv (v) (I-17) 
uy 


ae is given by 


a, fio ad oeua be seea ul) ih 


= ieee Ds [(pp]aq) + (q:pp)] + me Ras oe) 
P ap 


le 
(I-18) 


where (pp|qq) is the coulomb repulsion integral between 


{61-T) 
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atomic orbitals Xp and Xq! (q:pp) and (r:pp) are coulomb 
penetration integrals between Xp and neutral atom q 


anata lnus 
. =— = * = 
(q:pp) fos (11) ‘soul Xp (H) dv (I-19) 


In this method complicated atomic two center integrals 

were, calculated’ and evaluated explicitly. Only the integrals 
involving overlap densities between Seer neigh- 
bours were neglected. They applied their method to 

benzene. The method does not work very well. The singlet- 
triplet splittings. appear but are too big. The integrals 
are very difficult to compute. Moreover, the results 

of a calculation on the oxygen molecule using this method 


(7) 


are in remarkable discord with experiment. Tu Saeen 

order to develop a workable and useful 1n-electron 

theory tne, fol lowing, points, should be; considered: 

(1) mw-electrons should be treated separately from other 
electrons. 

(2) The orbital wave functions should be made up from 
m—electron orbitals on the individual atoms. 

(3) m-electron repulsions should be considered. 

(4) The treatment of nm-electron repulsions should be 
simplified. 


(5) The theory must have some provision to use atomic data. 
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(6) Empirical elements should be included as adjustable 
parameters. 
The Huckel method satisfies ep ance bo. and 
Goeppert-Mayver andeskiar method fails in 4,75, and 6. In 


(8,9) 


1953 a stheory WasPOULecOLward in, which ail) of these 


characteristics are present. The theory presented by 


(8) 


Pariser-Parr incorporates empirical elements in 


the antisymmetric molecular orbital configuration inter- 


action (ASMOCI) method. Pople‘) 


has introduced a set 
of simplifying approximations closely related to the 
Pariser-Parr approximation. Sometimes the Pariser-Parr 
method and the Pople method are combined as the Pariser- 
Hotere Loples(P-P—P) method. The main wleatures of «the 
Original theory are similar to those of the Goeppert- 
Mayer and Sklar method with the exception of the way 
some types of integrals are evaluated. The one-center 
core integrals are calculated in the same way as in the 
Goeppert-Mayer and Sklar method [equation (I-18)]. The 


int als §S and Bg are treated in the same wa 
integrals ae; Be y 


as inthe Huckel method, i.e. 


a (I-20) 
and 


core _)\a parameter when atoms p and q are bonded 
pq 0 otherwise | (eZ) 
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In order to reduce the number of electronic repulsion 
integrals the zero differential overlap (ZDO) approximae 


(10) 


tion was made, i.e. 


a) Xq (1) Ove (uj 0 form pid (I-22) 


This eliminates all electron repulsion integrals except 
those of the type (pp|qq) or (pp|pp) where the notation 


(pq/rs) means 


2 
= S 
(pq|rs) = fx a a x, (2) x5 (2) dv, dv. 


(IS es) 


10% 


Here the integrals are formulated in terms of real functions. 


In other words, in the ZDO approximation, the hybrid 


(pp|pq), exchange (pq|pq) and many-center (pq|rs) integrals, 


all vanish. 


The coulomb integrals, Yo! are computed theoretically 


or semi-empirically subject to the condition 


ae (I-24) 
ste R 
> > 
for large R (= [|r - eal When p = gq one can use the 
EO cn baw On Pariser ‘t) 
Paes Te ah (i225) 


where I is the ionization potential and A is the orbital 


GleEceron auLinity. 


(12) 


Mulliken has estimated Le@and A to. be” Lis22 
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ev and 0.69 ev respectively for carbon atom. Thus 


(pp|pp) = 10.53 ev. Using Slater orbitals with the 
orbital exponent Z = 3.18, an explicit theoretical 
formula (13) gives 

(pniipp) = 16. 93. ev (I-26) 


The applicability of the ZDO approximation can 
be justified by assuming that orthogonalized atomic 
orbitals (OAO) are used as the basis functions instead 
of overlapping atomic orbitals. 


The orthogonalization of a non-orthogonal complete 


set of functions can be done in many ways. One of them 
was proposed by Yara: and according to him OAOS 
can be represented by 
~1/2 
i= -yies / (I-27) 


where X and xy are row matrices and S is the square matrix 
of overlap integrals in the nonorthogonal basis. When 
the orbital basis is transformed according to equation 
(I-27) the matrix representation of any one-electron 
operator M will also be transformed according to the 


equation 


C23) 


A a he 


wl 





ayiT  .motn noes 102 ylevigooqaet va 3.0. 
edt dziw elsetidio oange pried ve £2,014 = (salen) 


{roijerosdd tioilgxe as ,8L.— = 3 Inanogxe tastes 


asvin **) otumscd! 


a 
ae 


=f) ve £2.8f = taqlaq) 


69 noiyemixorgqgs O08 are to vos itidaobteas old 
Simetcn fom isnopeds10. Jecie enrimrees vad Beilizeut ed 
bDeeteani eaenoitione} eieed sds #6 Bear szs (OAD) alasidse 
ne : 
Ast rdtb oreiter fn teqaizeve I6 : 

7 
etefonos [6nupofitgro~con & To notes erfagnopodsto sdiT 
moist Io sa .ayew yrem al enob sd aAsD anoitone? to gas 


yd betnossstget ed ass 


: 
2040 misfl oF pmibioovws fas ‘bwod yd beacqe2q. esw ; 





. , ‘ i 
(VS-z) a 2x A 
* 
xiiszem sisupe siz al 2 bas aepisdtem wot ote y Bae &£ etedw : 
: be ; 7 - 
qeqw .aised Lanopotitionon sd nt eferpacal geiteve Io _ 
nofjsups oJ paibroots batioteansz3 vi aiead invidso afd a 


noxrgosis~sno yrn io noissenaeeiges: xis¢em 6f39 {(TS-T} 


ent oF paibiroovs bamiciensis od cals [liw M z0Sez9ego . 





where ui refers to the matrix representation of the 
operator in the orthogonalized atomic orbital bas ts eae 
and matrices without superscripts refer to the matrix 
representation of the operator in the atomic orbital 
basis, {x}. It has been shown by Inga-Fischer- 

Hjalmars ‘19) Spot eranotommaciOng (l= 27)} tO (1-28). will 
be approximately equal to ZDO approximation. 

There are many choices which can be made for 
expandingemoOlecUlar Orbitals. These are (a) Slater-type 
orbital, (b) distorted atomic orbital, (c) combinations 
On) (a) "and™ (bp). 

In the Pariser-Parr-Pople method one assumes the 
second choice of the orbital where the explicit form of 
the distorted atomic orbitals does not need to be 
specified. Here all the parameters are used as 
empirical quantities except coulomb repulsion integrals. 
This procedure is justified because the operator involved 
is simple and the asymptotic form of the interaction is 
well understood. 

In the Pariser-Parr-Pople method both the ZDO 
approximation and the determination of certain integrals in 
an empirical way have been included. It has been shown 
that the ZDO approximation can only be applicable to 
orthonormalized bases. As a result the local character 


of various integrals will be destroyed. Hence the trans- 
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ference of numerical values from one molecule to another 
will be prohibited. However, Inga Pischer-Hjalmars ‘+>) 
has shown that the pertinent integrals are independent 

of the surroundings. 

Although the refined form of t-electron theory 
explains many of the properties of conjugated and aromatic 
systems it still cannot explain many chemical properties 
such as transition-state mechanisms in chemical reactions, 
steric hindrance and ring strain, etc. Moreover, the 
HMO and P-P-P methods were developed to study only 
planar conjugated molecules. Attempts have been made 
to use these methods to study non-planar molecules with 
go electrons but they are of limited success. An 
improved theory is needed in which calculations should 
be basedjon aljltvalence electrons. One of the first 


(16) who attacked 


successful attempts was made by Hoffmann 
the problem of treating all valence electrons ina 
quantitative manner with the application of the HMO 
method. In his method known as Extended Huckel Theory 
(EHT), all valence electrons are treated explicitly, 

all overlap @entegralssare calculated, only els electrons 

of carbon and heteroatoms are in the core and like the HMO 
method the Hamiltonian remains undefined. The secular 


equations and equations for orbital coefficients are 


given by the following equations: 
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lH - ES |= 0 (I-29) 
Eo wi . 
1%, a ES.) Cy 0 | (I-30) 


The summation is over all valence shell atomic orbitals. 


The matrix elements Hos are defined by 


Hy - fr, H e dv (I-31) 


The Hamiltonian H represents the energy of electrons 
in the field of a core made up of nuclei, inner shell 
electrons and the valence electrons. The diagonal ele- 


ments are determined semi-empirically and off-diagonal 


elements are estimated by dayoolniniedbaiesay Be! Forma. 
Thus 
Hae ee rT (I-32) 
UU u 
and 
Fpem eee (years can) S (ice) 
Lv 2 U Vie Cw 


WierteceKeis taken to be 1275. Off-diagonal elements 
Bicludesall PalLrs OL atomic orbitals. ~ Ther totadi, energy 


of the molecule is given by 
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Fuleisenotcerree [rom criticisms which are 
applicable to HMO theory. It does not predict the correct 
heat of atomization and in many cases the charge density 
Is far from reality. "The*fatvare"or-BAT led to the 
development of other sophisticated all-valence-electron 
methods in which t-o interaction has been taken into 
account. Among these methods CNDO (complete neglect of 
differential overlap) is the first theory proposed 


Dy AFOple et ayy?) (18) 


In 1967 Pople and coworkers 
also introduced the INDO (Intermediate neglect of 

differential overlap) method and this was followed by 
ee) 


Dewar and Klopman's partial neglect of differential 
overlap (PNDO) approach. In quick succession many other 
methods have been developed in this line in which the 
choice of parameters or method of parametrization is 
different from each other. As our work deals with all 


valence electrons, we shall discuss the methods in 
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Chapter II 


SCF THEORY AND ITS DEVELOPMENT FOR 


T-MOLECULAR SYSTEM 


In order to discuss the all-valence-electron 
semi-empirical SCF theory, it appears most profitable 


to first outline the formulation of the SCF MO method 
(20) 


(9) 


by Roothaan and its application by Pariser-Parr ‘8) 


and by Pople using various approximations which are 


in the area of semi-empirical choices for the integrals. 


Hartree-Fock-Roothaan Method: 


The central problem in the elucidation of the elec- 
tronic structure of N electron systems in a stationary state 
by Poven teehani cal means is the solution of Schrodinger 
equation 


HW, = EW. (II-1) 


where H is the total Hamiltonian operator of the system, 
WY is the wave function and E the energy eigenvalue. The 
wave function depends on the electronic coordinates each 
consisting of a space coordinate and a spin coordinate. 


It is well known that the solution of equation (II-1) 
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involves considerable difficulties. In order to obtain 
a good approximate solution of equation (II-1) the 
Variation principle is to be used, Expressing the total 
wave function for N electrons as an antisymmetrized 
product (AP) 


v,)) v, (1) soe C1) 
; (II-2) 


For closed shell systems the molecular spin orbitals 


(MSO) are given by 


Wohi oe Mae Ue ma ea 


With 2n electrons the antisymmetrized product may then be 


Written’ as 


i 
Yo = [(2n) 1] 2), (1)0(1)$, (2) 8(2) ... 


wets >, (2@n-1) a (2n-1) o,, (2n) B(2n) (11-4) 


The MSO's are assumed to be orthonormal, 
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An important implication of the above equation is that 


the AP given by equation (II-2) is normalized i.e. 


[wy dr =1 Chie) 


The electronic energy for an electronic state represented 


by the normalized wave function is given by 


E = fFnvar (II-7) 


where the total Hamiltonian, defined earlier by equation 


(E—- 13) fpebsn tas. Ws 


H HEV TLV 
Pes 
with UG) LA eee ae (II-8) 
Z a Fau 


Substituting the wave function equation (II-4) into the 
energy expression (II-7) and performing the necessary 


algebra, we find the energy of a closed shell AP to be 


E = 25H, + 24 (255 .-K, .) Giro) 
i ay 


where the nuclear field orbital energies Hay coulomb 


integrals ay and exchange integrals Shs are defined by 
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Hone = fF, H (u) 6, (av, (II-10) 


me hy Seimei sere si me ae a 


cae 
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1 ) (II-11) 
Ks aR eR soffF, °. (v) (ex >, (u) o ONES ete 


eromvequation (11-11) it as clear that 
a eh (II-12) 


It is useful to define a coulomb operator J; and an 


exchange Operator K; by the following equations 


| $.(v) o, (v) 
3,6) = |S o(u) 


Diy 


(II-13) 


These operators are hermitian and linear. J. is the 
operator expressing the potential due to an electron 


Gistiabutedvin space with the density Peal However K, 


has no classical analog. 
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In order to obtain the AP for which the energy 
of the system reaches a minimum, the expression (II-9) 
is to be minimized by varying the MO's within the limits 


that they form an orthonormal set such that 
a = 8), | (rr) 


When each MO os varied by an infinitesimal amount ob 

the energy is also varied by an amount SE. The 

condition for E to reach a minimum is that SE = O for any 
infinitesimal eon of the o,'S. This econd1c1on) is 


given by 





Heer OA(Z0NSK (Gye = ee be 
} ae BSS ote 
J J 
(II-15) 
Hy 43 8 20K O KG PE eree ey ° oe 
( . ifs o. ety 
3 J 
where ae avelagrangian multipliers. (lt canebe shown 


that Lagrangian multipliers are the elements of a 
hermitian matrix. As a result the two equations (II-15) 
are equivalent. The equations are known as Hartree-Fock 
equations proposed simultaneously and independently by 


(21) (22) 


Fock and Slater 












vergar 22) rinidw 1o2 TA) ang nieido ot Tebpe ae i 


¢ & i 

iGeET) «A rove. Ody  sttin kd 5 estivess Seg ea Oe 
a, 

atimii ori: coe ‘i ai pilyaey yd Cesantnd ga Oe Bae 


ue $45 j 1 ’ so (Gl mo voaly Send2 
a. 
P e 


im | x8 | eice .@ OM doses ca ee 
a 
. P, : ) — 
ho 46! 218 24 YRaeae Pi oe - 


wns xo) 7 “ | Jeo, G? 3 40% aoldieie ee 
a 
- ' Jvtsivev tenseed nigel 


¥ gd sovip. - 


- 


TT, 
“ 
- 
i 
Cc 
¥ 
Sail 


ett: 
A! 
— 
ts 
' 
he 
vw 
~~ 
ian 
+ 
~ 
———— 


nrwonie SuttiGo . Sergi Ley ALLE] GIB seo 


; cf . 
& 36, ci jcvnet> of). Sue Bis iig-y lun nae 


| GELAOE) Grolsaups owt of? -tivesn & 2A eee 










Bo Myosd/ eae, sive Gare oT 


21, 


We now define the Hartree-Fock operator 


Revs Horr (20 .<K.) (II-16) 
chee 


For the best MO's to satisfy the equation (II-15) we 
can write 
Fo; = 2 ou (ed 
eee) e- 
It is known that without changing the physical content 


of the total wave function (II-2) we may cast the equation 


(lisi7) into the following form 
(II-18) 


The set of equations (II-18) is most commonly known as 
Hartree-Fock equation. 

The general procedure for solution of the Hartree- 
Fock equations is to choose a set of trial solutions, 
from which coulomb and exchange operators are formed and 
the Hartree-Fock Hamiltonian operator is calculated. Then 
the equation (II-18) is solved for n lowest eigenvalues 
and the resulting o,'s are compared with the assumed ones 
and a new set of o,'s is chosen from the knowledge of this 
comparison. This process is then repeated until the 


energy value of the calculated and assumed o,'s agree, 
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This method for solving the Hartree-Fock equation is 
known as the Hartree-Fock self consistent field (SCF) 
method, 

Although the Hartree-Fock method gives a useful 
approximate solution to the Schrédinger equation, it has 
been used in its original form of integro-differential 
equations almost esaailetnckiee forvatoms., For molecular 
systems which lack spherical symmetry, it is almost 
impossible to solve the equation SrLectively In order 


ReO,Overcome this ditficul ty Roothaan ‘7°? 


has developed 
a method to obtain solutions to Hartree-Fock equations 
avoiding numerical integrations. 

Kee {x3 represent a complete basis set of 
functions which are normalized although they may not be 


orthogonal in general. Each Hartree-Fock orbital can 


bewexpanded in terms of their basis set 


5 = 2 Sui = xX C. eee @ © @ (II-19) 


where C; is the column matrix of m rows such as 


C;= and xX = (Xy Xo 22 +X,) 
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1a {x3 constitute a complete basis set then fo, 
are given exactly. However in practice the expansion is 
truncated to m finite members of the basis set. The 
condition for constructing n linearly independent 
solutions is that m 2n. 

Ite {>} form an orthogonal set then it requires 


that 
>» ore S ie ae (II-20) 


where by is the Kronecker delta and ane is the overlap 


integral for atomic functions defined by 


Hi] 
= 
J 

Be 
(ES 
O; 
i 


Shae 
1) 


<a De Cee =) Oo! (II-21) 


where 


and S is a hermitian matrix the elements of which are 
overlap integrals represented by the equation (II-21) and 


diagonal elements are unity. 
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24. 


Let us define the following matrix elements and 


matrices 
= Cx od. 
iS elias ies (II-22) 
* * 
(II-23) 
* * 


In order to obtain the best LCAOMO's a variational 
treatment is to be carried out. When the vector Cc. is 
varied by an infinitesimal amount 6C., the coefficients 
Che are also varied by a small amount ei The condition 
SOEpthneGpenergy steosbesarmingmum is GEe=203 aThe Euler 


equation for this variational treatment assumes the 


following form 
2 Ges = Ss C. oa (ides 2) 


or in a more compact notation. 


RB €y- 15aGa¢ (II-25) 
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The equation is known as Roothaan's equation. If ¢€ 
represents a diagonal matrix with real diagonal elements 


5 the equation (II-25) takes the form 


Fe, = ene OMG. (II-26) 


The non-trivial solutions of equation (II-25) are obtained 


by solving for m roots of the secular determinant 
dete ers | =O (II-27) 


Since the matrix F depends on tc, } the secular equations 
are non-linear and should be solved by iterative 
techniques. This can be done as follows. As a first 


va) 


approximation some matrix € is assumed and 


Weis ze0 sco. corstruct pll), then equation (II-27) is 


solved to obtain a first improved matrix ¢ '2) , uy om} 
process is continued until the matrices between two 
successive iterations agree within a specified limit 
when the solutions are said to be self-consistent with 
respect to elements of matrix €. Since the basis 
orbitals {x} are not varied this method does not give 
true Hartree-Fock solutions unless the basis used is 


large. The lowest n roots correspond to those occupied 


by 2n electrons describing the ground state of the atom. 
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The rest of the solutions are known as virtual orbitals 
and may be used to construct excited state configurations. 
The energy can be improved by the configuration inter- 
action method. 

Ab initio wWartres~Fock=Roothaan molecular 
calculations are limited at present to moderately large 
molecular systems. * The difficulty is inthe evaluation 
of the integrals of the type (sya ee) because the number 
of such integrals is very large and increases roughly as 
the fourth power of the number of functions of the basis 
set. Moreover, the determination of these integrals 
presents technical difficulties for functions which are 
Orbitalsuor three, or four ditferent atoms... In addition, 
it is seen that the results of such calculations are 
frequently not in good agreement with experiment. Hence 
it seems that it is desirable to have methods which give 
better agreement with experiment and which are simple 
enough to permit extension of our calculations to larger 


molecules. 


Pariser-Parr-Pople Method: 


The essence of the Pariser-Parr-Pople (P-P-P) 
method for obtaining the wave functions and energies of a 


given 7 system is as follows: 
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(a) The Hamiltonian operator for the 7 electrons of a 
given molecule may be represented by equation (I-13). 
(b) The molecular orbitals may ba eaohase aa as linear 
combinations of atomic orbitals as in equation (I-5). 

The energy of the system containing 1~- electrons 
can be determined by using the Hartree-Fock-Roothaan SCF 
method. For this purpose we need to know the value of the 
electron interaction integrals. The applUcel 1 Onior tiie 
ZDO assumption reduces the number of integrals from the 
order of n* to ne where n is the number of t-electrons 


of the system. The matrix elements are TepLresented by 


the formula 


Sg ss 1 Pu Von e 2 Ry Ge (II-28) 
2 hata s cae aN (u+v) (11-29) 
where o,, = (x, (0 Bieay Ge hg, (2) (II-30) 
Bay Be Six yloghe Beet js ane) (Tt-se) 


ae is the bond-order matrix and ins is the coulomb 


repulsion integral. 
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Now OF is the energy (kinetic and potential) of an 


electron in the orbital Xu in the framework field. If 


we write the potential energy term as VEV tVo+ee eee We 
have 
gee 
= -= VY + 
hy Gale h its Myles Og yy bg 
= + II=32 
U, DE the, lwglx ad ( ) 
V+U 
apds ifs potelxela ——sav2sve tv. clog it ie ala ale ei a3) 
uv ul ee Thy ee Lee pu y 
PFH,YV 
Thus 
ab 
Ree = yey +> -P + - 
si Sa : es Fe, Oly Py? (II-34) 


Uys which is regarded as an empirical parameter, can be 
taken as constant for a given atom. If Xu is an eigen- 
function of the isolated atom then oF is the energy 

of an electron in the orbital Xi and this may be 

taken to be equal to the valence state ionization 
potential “7 of the atom i.e. the ionization potential of 
an atom in its proper hybridised state in the molecule. 
For a given pair of atoms and given bond length the first 
term in Bay can be taken as constant. Since the potential 


due to a distant nucleus having unit charge is equal and 
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Opposite in sign to the potential due to a distant 


electron density pay BterOllows that 


OG Pele -Z,, (uv | pp) (TI-35) 
where os is the charge on the core at atom P. 
According to the zero Overlap approximation two 
electron integrals of the type in equation (II-35) are 
equal to zero. Hence the sum over distant cores in Bay 
can be neglected. Thus By can be taken as an empirical 
parameter and is constant in a given type of bond; 
By is assumed to De zero for non-bonded atoms as in 
Huckel theory. But it should be noted that the 8 value 
used in Huckel theory for a particular type of bond is 
“not the same as that needed in P-P-P method because in 
Huckel theory 8 compensates the neglect of all electron 
repulsion terms. 


Adopting= the approximation. (11-35) for'the* case 


H=v, we have 
Caos bee Zoay (II-36) 
H HL ne, p Up 


If we use the value of Os in equation (II-36), Eat takes 


the form 
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In the ZDO approximation Pople used expression (II-29) and 


(II-37) for the construction of SCF molecular orbitals. 


(8) 


The ZDO assumptions suggested by Pariser and Parr 


and by Pople”? may be summarized below: 


(ES) oN (II-38) 
(EB coal Oh ea (II-39) 
A ere) + 0 (where y and v are 

neighbors) (IIT-40) 

(u[Hoo+e!¥) = 0 when y and v are non- 
neighbors (II-41) 
= 6 II-42 
(uv | kn) oy ean ee ( ) 


where 


2 i 2 
ia = fix Fp LX, (2) | dv, dv, 


The ZDO relation of equation (II-42) was suggested by 


Parr previous to the assumption (II-38) to (II-41). 


30 














Se gg am we Dyini at? Fy oe ™ 
| | 7” | ee ena 
| we). 
fe (@$-I1) noiadetyxe been siqot noltsmixotags OGk ad? ab 7 
-disdidso xalivoslom Y2 to nolssuxtenco edd we (V&~IT) 


9) sad bis iseizei yd besesepouse enorsqmuses ods’ oft : 
twofed bes bisa od nal) oltgos ya bos : 
(oe~02) | ii = (vise 
(Q&=11) | > + Tae) Tis 
S38 v brs yy exzerw) 0 > («1 Seegtthae ; 1 


(Ob—ZL) (axddidntons 


“Men exe v bas yu vodwd = 
‘ 44 (v1.5 gabe: 


(fa-I1L) axodripisa . 


(Sb-ZT) 





Sly 


These equations (II-34) to (II-41) were introduced by 
(8) (9) 


Pariser and Parr and Pople 


(15) 


as useful simplifications. 
But Fischer-Hjalmars has shown that these approxima- 
tions were largely justified if the basic atomic orbitals 
are taken to be Léwdin orthogonalized atomic orbital 

(OAOC)* instead of Slater erbitals, 

In t-electron theory the core resonance integral, 8, 
is treated as an empirical parameter and is determined to 
fit an experimental property. For example in recent 
yeates there are two groups Of @ values for a C-C bond (1.40 
Dj oundeinethesliterature. ~One, higher value —2.39 a 
is used for the prediction of spectral transitions and the 


other lower value (-1.75 ev) (23) 


utilized to predict the 
ground state properties such as heat of atomization, 
resonance energy and bond energy, etc. 


There are several methods for evaluation of 


resonance integrals 


(a) To fit the lowest singlet~-singlet transition 
energy. 
(b) ChungeandaDewar | sanduno sandewnitehesd. a) 


evaluated the integral using a hypothetical thermocycle. 


(c) In the 8 variable method as proposed by Yamaguchi et 


Fr cao two center core integrals Bye the bond length 


Gs and the two-center repulsion integrals Nee between 





neighboring atoms are adjusted at each iteration 
according to the following equations until the self- 


consistency was achieved. 


° | — 
Diy (Ae) ey. 520 0.186P, , 
(II-43) 
cae = 8 exp | a (Zo-F yy) | 
Ly is the C-C distance for benzene, a=1.7 Rea and By is 


the resonance integral for C-C bond in benzene. 
(d) Resonance integrals may be related to overlap 


integrals by the equation 
8 = -kS (II-44) 


where S is the overlap integral evaluated using Slater- 
zener orbitals (z=3.18) and k is a proportionality 
constant; for carbon-carbon the value is taken to be 
eyon7 os 


Two center two electron integrals (pyul|vv) aces 


calculated in the following ways: 


(8) 


(a) By the charged sphere approximation ‘ 
(b) Pople method: 
(uu}vy) = BEB OOTE ta sy Sia ee 
“ab (II-45) 
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(Cc) Method due to Mataga and Neehimotoreu 
CO = = oeso where r is in A° 
ie (Tr 169 
ik : : 
= where r is in a.u. 
a+ x 
uv 
where : a 
ee 14.3968 ey He ekg 
Ll Gin te il Li 
5 : A) + 5 (I, A) 
(II-47) 
| : 
e Tina s ls 
1 Rs 
pea oo etn (I-A) 
(da) Ohno's method '2°) . 
(uul[vv) = etl 1/72 where r in A? 
2 2 
reels ) (II-48) 
iit : 
= a where r in a.u 
ie. 2 
(ar Os) 
UV 


where a is defined as in equation (II-47). There are 


many other formulas for the evaluation of y integrals 


(le) (28) 


but Mataga and Nishimoto's method and Ohno's method 
are most commonly used. It is seen that for the study of 


ground state properties of molecules Ohno's method is 
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better than Mataga-Nishimoto's method whereas Mataga- 
Nishimoto's formula is superior in the context OFeC TL 
studies, 

It would be an interesting invention if we could 
construct a simple formula which would reproduce both 
Ohno's and Mataga-Nishimoto's values and also Roothaan's 
theoretical values of coulomb integrals as special cases. 
In the following we introduce such a formula which will 


be discussed in more detail in Chapters Ly... 
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Chapter III 


ALL-VALENCE-ELECTRON METHODS 


The importance of molecular orbital calculations 
fata be much increased if the methods which have been 
described in Chapter II were extended to all valence 
electrons instead of being restricted to the level of 
m-electrons only. Such an advance would permit not only 
a full treatment of o and ot electrons in planar molecules 
but would also allow the calculation of a large number 
of molecules where o-1T separation is not possible. 

Since EHT is not free from disadvantages similar 
to the Htickel theory, a theory concerning all valence 
electrons is needed in which electron repulsions have 
been taken into account. In the following we would like 
to discuss some of the methods in which all valence 


electrons are included. 


The CNDO Method: 


The CNDO (complete neglect of differential overlap) 
Ce) 


method proposed by Pople et al assumes that the 


peocuce ofsorbitals XX is always zero provided that 
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Ke a ciatal X,, are ditferenty! According to this approximation 


u 
the overlap integrals SH between any pair of different 
orbitals is zero except in shies (II-6) and all integrals 
(uv|Ac) are zero unless u=v and \=c0. The integrals 

(up fAdA) are written as vain and are assumed to depend on 
the atoms A and B to which Xa and X belong “And Hot on 

the type Of orbitals. Thus all the various ei between 


atomS A and B are approximated by I The rap are 


AB* 
approximated in the original CNDO calculation by the 
purely theoretical values of the integral (8.8, |S,8,) 
Or ethe valence shell s orbitals. 


With these approximations the Hartree-Fock matrix 


elements become 


aN ios aay ‘3 (Pan FP iy) “798 f STA. BB AB AB) (TIT-1) 
— a — 
aK = gy x uy PAR (uv) (ITI-2) 


where the atomic orbitals My is connected on atom A and 
X On atom B. In deriving these equations Pople separated 


the core matrix elements Sen Thus 
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(III-3) 


The diagonal elements Bey Maybe split up anto two parts: 
the first term aes is that part of the diagonal matrix 
elements involving the one-electron Hamiltonian 


containing only the core of its own atom and terms, Vapi 


containing the interaction of the electron in X, On atom 


A with the cores of other atoms B and is approximated by 
eae 
the relation fee . dv. The valuet of U is 
ae HU 





approximated from atomic spectroscopy. The Pos, are the 


components of the charge density and bond order matrix 


P 2 doc G (III~4) 


and ot the total charge density on atom A is given by 
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Finally the off-diagonal elements Abe are taken as zero if 


Xu 


different atoms A and B, then 


and x, are on the same atom. But if Xerana Y "are ou 
Vv pees Vv 


wed Dake: ° o 
H St (Be +8e ) “bes | (III-6) 
where Bx and BS are treated as adjustable parameters. 
Under the CNDO approximation the total energy of the 


molecule can be written as a sum of one- and two-atom 


terms, 
Total Hee — AB (III~-7) 
where 
A A A 
Ss é ih hice) 
Paes aL, eter wu a De os ee oo) ay) aA (III-8) 
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The interaction energy of the cores of atoms A and B is 
~l 
approximated by Z,2Z.R where Z, and Z,, are the charges 


AB AB A B 
of the cores in units of +e, 

One of the principal defects of CNDO/1 is that 
for diatomic molecules the calculated bond lengths are 
too short and binding energies are too large. This is 
due to the 'penetration' effect of electrons i.e. 
electrons in an orbital on one atom penetrate the shell 


of another giving rise to net attraction. Writing the 


equation (III-l1) in the form 


eye) (III-10) 


the last term may be taken to be the contribution of 
penetration integrals to Pais In order to eliminate the 
above mentioned difficulties Pople and Segal‘? ina 
later paper modified the CNDO method known as CNDO/2. In 


this method the penetration terms in equation (III-9) are 


neglected by putting 


V set ea (III-11) 
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Although there is no complete Ousti fication for this, 
Pople and Segal ‘17) have suggested that Meg lectsof 
Overlap causes errors which are Opposite in sign to the 
errors due to the neglect of penetration, Thus equation 
(III-11) corrects the CNDO/1 procedure in the raghnt 
direction. 

Another modification made in CNDO/2 is the method 
of estimating the local core matrix element Eee In-*CnDpo/ 1 
this was done by the relation 


Sal =U + (2,-1)T, (III-12) 


ei HE A 


where a is the ionization potential in the appropriate 
average atomic state and orbital iy belongs to atom A. 
An alternative procedure is to use atomic electron 


affinities, 


ie ly = U teen eats Ls (III-13) 


In the CNDO/2 method the integral ha is taken as the 
average of the ionization potential (I) and electron 
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Thus 


(III-14) 


and 


(Cen) 


Using equation (III-11) and Get r=15) the basic equations 
for the Hartree-Fock matrix in the CNDO/2 formalism may 


be written as 


bras | es wi 
Se a, (ai +A) + an Zn) = ( te Dan + 
py) (Cee eG (ITI-16) 
Bia BBB) AB 
es 1 ieee ut i: 
a x S) (B2+8s) ae aitine SNS (III-17) 
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Invariance Problem: 


Pople and co-workers ‘17) pointed out that in 

a full Roothaan-SCF-MO treatment the two-—center integrals 

are independent of the choice of axis but this is not 

so in simplified versions because these integrals are 

dependent on the choice of coordinate system and 

hybridization of the orbitals. Hence in the zero 

differential overlap methods the results of these 

integrals ape required to be invariant, to two types of 

transformations: 

ci) _ Coordinate Invariance - invariance under rotation 
OF the lcoordinate axis. 

C2.) Hybridization Invariance - invariance under 
rotation of the hybrid orbitals. 


We shall discuss this invariance problem in more detail 


mnechaptercy. 
Evaluation of Integrals: 


In the CNDO/1 method the overlap integrals are 
calculated theoretically. The electron repulsion 
integrals PAB are calculated as the two-center coulomb 
integral involving valence s eihermere which have also 


been evaluated theoretically. 
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Defects of CNDO Level Methods: 


Although the CNDO methods proposed by Pople et 


, (47) 


have been used successfully to calculate the 
properties of many species, they are not free from 
Shortcomings. One of the serious defects of the paramet- 
rization of the theory is that the total energies of the 
occupied orbitals are uniformly too negative and the 
energies of the virtual orbitals are similarly displaced 
downwards. As a result the ionization potential and 
electron affinity are too large by several electron volts 
for all species REO Moreover these methods do not 
Givesthe correct energy,or correct bond length, relation 
although they are good for bond angles and charge 
densities ‘1”) , The heats of formation calculated by 

CNDO in its original form were not in agreement with 
experimental results. As a result several authors tried 
to reparametrize CNDO methods in order to correlate the 
calculated and observed properties. untsucpea tiie by 
modifying Pople's parameters, was able to calculate 
structure (bond lengths and bond angles) and heats of 
atomization. He reduced the 8 proportionality constant 


and changed the ionization potential of the carbon atom. 


He also used a linear relationship between calculated 


43. 


ieee , | rm. yor 
LiTaa) Stusou 


‘ i 
sit bsovbs1 oN .nozzesio 


‘Is strage me 
es 





_ 44, 


and experimental heats of atomization i.e. in order to 

get good agreement between calculated energy and observed 
energy of molecules he multiplied the theoretical value 

of energy by a scale factor. It is not clear how 
reasonable this procedure is. Although Davidson et aah 
used Wiberg's parameters in their CNDO calculation they 
did not use his linear adjustment of energy values; 
instead they used the calculated energies which are 
obtained from the usual CNDO procedure. Fischer and 


Kollmar (31) 


reparametrized the CNDO method and obtained 
good heats of atomization, equilibrium geometries and 
force constants. For penetration they have used the 


formula 


= a [Ser Re Gik7; te = 
VB = Zn (1-a) ae FOvR + L/ ix (ITI-20) 


where a=0.22 and Un is the orbital exponent. Although 
the equation (III-20) works well when R is small it is 
doubtful how far this equation will be successful when R 
is large as in large molecules. Theoretically, the 
equation (III-20) is not sound because when R + » V ae 


AB 


o whereas in reality when R> % VaB oO eels teruecO 


obtain good results, they have changed the functional 


forms for resonance integrals and certain nuclear 
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attraction integrals from those in CNDO/2. They have 
adjusted a total of eight parameters and the requirement 
of hybridization invariance was relaxed. 


(19) and Dewar et ai t9a) 


Dewar and Klopman 
Calculaceds heats Of formation of a large number of 
hydrocarbons by using PNDO (Partial Neglect of 
Differential Overlap) and MINDO (Modified Intermediate 
Neglect of Differential Overlap) methods byeaintroducing 
rather artificial nuclear-nuclear energy repulsion terms. 
Although the results were in good agreement with 
experiment, how such an ad hoc approximation is 
justified in order to get good heats of formation is far 
from obvious. It is also seen that such an approximation 
yields unrealistically short values for bond length. 
HPAUSeUnbe MiG resultwis tO Sacrifice One at the cost of the 
other. 

The choice of parameters or mode of evaluation of 
certain integrals discussed so far is useful only for the 
study of ground state properties of molecules and cannot 
DesuLuta 2Z6dsineg predicting theywspectral transitions. 

Thus the conflict between ground state properties and 
spectroscopic properties also exists in all-valence- 
electron methods as well as in the MO theory for m7 systems. 


Thus it would not be good to transfer certain parameters 
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for the study of certain properties to the Stuay or 
others because they account for different physical 
properties. 

Del Bene and saree 32) uSing limited configuration 
interaction studied the spectra of benzene, pyridine and 
some related molecules, They also used different 8's 
for o and tm. Although they are successful in predicting 
the spectra of some conjugated systems, the method fails 
to predict correct heats of formation and molecular 
geometries. They adopted the CNDO method for the 
Calcudation of the energy levels of those molecules and 
their aim was to gain information on the effect of the o 
electrons on the tm energy levels. 

So far in all these methods a minimum basis set 
of valence shell A.O. has been used. Hence these methods 
could not account for transitions involving a change in 


(33) 


principal quantum number. Salahub and Sandorfy used 


higher orbitals of H (Ho Si pte Wie Gahash Ae Seep dae viel, Gee ele 


2p 
CNDO method for some paraffins and obtained reasonable 
interpretations of the main characteristics of the 
observed electronic spectra. 

Hence we can say that CNDO calculations which give 
good results in molecular geometry, do not give good 
results in heats of formation and spectra; again when it 


gives good results in spectra it does not predict ground 


state properties correctly. 
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Chapter IV 


KLONDIKE FORMULA 


There are several different formulas commonly 
used for the computation of coulomb integrals in semi- 
empirical molecular orbital calculations. 

As described earlier, people have more or less 


settled for Ohno's recone ee 


(27) 


or Mataga-Nishimoto's 
formula in the case of the t-electron approximation 
(p-p-P theory). In the development of the all-valence- 
electron approximations the situation has been far from 
settled; some prefer using Roothaan's explicit theoretical 


(13) 


formula , while others carry over their habit acquired 


in the t-electron method to the new stage and use either 


Gnnots \-°) 


or Mataga-Nishimoto's ‘?7) formula. Now that 

it has become clear that the choice of the formula for 

the computation of coulomb integrals constitutes an 
important part of the semi-empirical adjustment of the 
individual method, it would be convenient if we could have 
a suitably parametrized formula which would cover conti- 


(13) (28) 


nuously the results of Roothaan's Ohno's yp and 


(27) formulas. “We have succeeded in 


Mataga-Nishimoto's 
devising such a new artificial formula for coulomb integrals 
in our present work. This we call the Klondike formula, 


which has the following form, 
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GIRLY = nepali (IV-1) 


where R is RAG the distance between two centers yu and \v, 


and 


where At and K\, are the klondike parameters for yu and vy 
centers respectively. In the present work we have 
restricted our considerations to the case Sie ia 

When k=0, the above coulomb integral formula equation 


(IV-1) is identical to Mataga-Nishimoto's formu! a 


The value of oy can be determined in the 


following way: 


when R = 0 


(uufvy) = 


when the two functions are identical i.e. when y=v 
= 2a 
(un | uu) i 


For Roothaan's method the value of (uu | up) is determined 
explicitly. For Mataga-Nishimoto's method or Ohno's 
method (uu|upy) is equal to I-A. 


Thus 
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a= 5 (TA) = Gul) 
We are designating the method using the Fences formula 
as CNDOK and the method using Roothaan's formula as 
CNDOR, 

In our parametrization we have found that when 
K=0.4, the Klondike formula corresponds to Ohno's formula 
and when k= 0.67 it corresponds TO yRoothaan's formula. 
and 


Table (1) compares the values of T r 


bsilsi~. 282s¢ 
iciae calculated by the Klondike formula and Ohno's 
formula. Figures 1, 3, 5 represent the variation of [ 
values with R. From these figures it is seen that in all 
cases the curve of the Klondike formula agrees well with 
tiwat-of Ohno; In the case’ of eA deviation between the 
Klondike and Ohnocurves is appreciable at small inter- 
nuclear distances. This is not significant because 
molecules do not exist with such small internuclear 
distances at equilibrium. 


Lab bews(2Z)a-cOntains the values. of 1 i 


ie Sie oo ss 

and losis calculated by the Klondike formula and Roothaan's 
formula and figures 2, 4, 6 are also drawn for comparison. 
It is seen that for ee the two methods agree well at 
abl internuclearpdistances, ~Buteron Toa? Po g28! large 
deviations are observed at small values of R. As stated 
earlier this has little effect in the calculations because 


no molecule exists with such small internuclear distances. 
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In the main body of the present work, Roothaan's 
formulas are used in accordance with the widely accepted 
CNDO/2 method. However, we shall also describe our 
preliminary attempt to use the Klondike formula in 


conjunction with the CNDO approximation. 
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Chapter Y 


A NEW PARAMETRIZATION OF ALL-VALENCE-ELECTRON THEORY 


In all-valence-electron methods at the level of 
CNDO approximation, the parametrization and choice of 
parameters are of central importance. A particular 
choice of parameters can predict certain properties of 
a molecule which are in good agreement with experimental 
results. However, use of the same parameters may cause 
other: predicted properties to disagree with experimental 
values. One of the most serious examples is the dilemma 
between good bond lengths and good heats of formation as 
mentioned in Chapter III. In this work we have made an 
attempt to parametrize the theory in a way to obtain 
good correlation between energies and bond lengths of 
some hydrocarbons such as methane, ethane, propane, 


butane, ethylene, allene, acetylene and benzene. 


Description of the Method: 


The present method is similar to the CNDO/2 method 
which has been described earlier. Here we are going to 
discuss the essential features of the present method. 

In the SCF-LCAO-MO theory the molecular orbitals 


are determined by solving the Hartree~Fock-Roothaan SCF 
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equation 
FC = SC€E (V-1) 


where € is the matrix of the coefficients Oe ae E is the 
eigenvalue matrix of orbital energies and S is the overlap 
integral matrix, 

The elements of the Fock matrix Sey may be split 


Hpeincormtwo parts 


ap = aa + Sa (V-2) 
where 
A |: wiy2_ _ 
Sey =fucay 5V 2, WROD) Woks, (V-3) 
(one-electron part) 
and 
ye \ if 
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The diagonal elements rd may be divided into two parts. 
One is an atomic term and the other is the interaction 


of the electrons in X,, with the cores of all other atoms. 


| rod | 
ot ie Ta Ae SOE (V-6) 


The off-diagonal elements Se are given by 


= a ° ° 
H 7 leven Say (V7) 


We have approximated the second term in Aa by the 


relation 


Gi pmmniitaez (i-Xe yaad (V-8) 


where Zn is the core charge on atom B equal to the 


nuclear charge minus the number of inner shell electrons, 
PAB is the electron repulsion integral and i and w are 
global parameters which we introduce in the present work. 


The justification of the formula in (V~8) will be 


discussed shortly. Since 


Say = tn (except in te. 
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where the summation is over all basis functions and 


ee | 
eae’) (V~10) 


Here Sigh refers to the diagonal elements of the bond 
order matrix i.e. the electron density in a given 
Orbaca li, 

In order to justify the modifications in the CNDO 
method introduced in this work and to evaluate certain 
parameters and integrals we would like to make the 
following comments: 

Ge) In principle, the one-electron part can be 
computed within a reasonable time, with the possible 


exception of three-center integral 
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If uw and v are on the same atom and ul or Vv represents one 
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epson ee | Mui eyes 
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C+A,B 


(iBone A, ©) on 7B) 


If appropriate expressions are given to iv,t, Au and ames 


may be computed explicitly, but we introduce a series of 


approximations. A primary reason for this is that we 


attempt necessary adjustments to the drastic approximations 


in the two-electron parts. Essentially, the one-electron 
part of the total energy is negative and the two-electron 
part is necessarily positive because they correspond to 
the repulsive electron-electron interaction energy. 

Thus any change in this positive contribution must 
definitely be counterbalanced by the negative 
contribution from the one-electron part. In order to 
incorporate this feature into the formulas we have 
introduced the approximation 
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Where in the present investigation \ and w are treated 
as global parameters. While this explicit formula is 
something to be debated, the logic behind it as 
described above is straightforward and seems to be 
difficult to reject. Another obvious compensation we 


may make, would be to drop 


7) (u/Vo] v) 


C+A,B 


in He because the corresponding three-center two-electron 
integral, say (HaV_ lo Qo) is also dropped by the ZDO 
approximation. 

(2) In the P-P-P method, various approximations have 
been used for (yuluy) and (uwul/AA). It may be said that 
reduced values for (un | un) arise from the empirical 
evaluation and for (uu | AA) convenient approximate formulas 


are used, the most popular ones being Mataga-Nishimoto's ‘#7) 


and Ohno's 28) | 

In the all-valence-electron method, the situation 
is somewhat different. Pople and his school make use of 
explicit Roothaan formulas for coulomb integrals. Clark 


(34) 


and Ragle contend that one should not use Mataga- 
Nishimoto type approximations because they would intro- 


-duce some inconsistency in the CNDO method. However, 
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other workers such as Del Bene and gafre 32) Dewar and 


Klopman ‘9) 


use those approximate formulas for coulomb 
antegqrals rather freely. | 

In the present work we have attempted to LALty 
these two distinct approaches by introducing a new 
approximate coulomb integral formula which we call the 
Klondike formula, as aus, BRR earlier. 
(3) Yee are the atomic matrix elements of the one- 
electron Hamiltonian. It includes the core potential of 
the atom to which Xi Detongs. It can be calculated from 
approximate atomic orbitals or it can be estimated semi- 
empirically from experimental data on atomic energy 
levels. It can also be used as an empirical parameter. 


Pople ‘17) 


approximated it from atomic spectroscopy. 
According to the differential Overlap approximation all 
Ereceeoulc states of a given configuration have the same 
energy, thus 


E (X, 2s™ap") =mu +5 (mtn) (mtn-1)R, +C  (V-14) 
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+nU 
Ss 


Z 8 


where term C comes from the inner electrons. In order 


toucalciuilate Sn from the experimental data, it is 
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therefore required to consider an average energy of all 
Wicestaren eC Petivesconriguivatton. “Lnwonderstogcalculate 
the single Ney the proper ionization process should be 


considered, 


m—1 


Io (X, 2s"2p”") E wea 2s 2p") - E (X, Demons) 


= -U - (mtn-1) T (V-15) 


An atomic orbital has the tendency to acquire and 
lose electrons and this tendency of the atomic orbital 
can be accounted for by the molecular orbital theory. 


As a result in CNDO/2 the following equations 


Ba = on + Zy i (V-17) 


have been averaged to calculate Sie Thus 


- 2 # = = 
re AN Rete ON aD iy ic) 


It is quite clear that the core parameters of CNDO/2 are 
uniformly too negative whereas the parameters of CNDO/1 
are less negative than CNDO/2. We have used oon as an 


empirical parameter in our calculations. 
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(4) If we can write molecular orbital v, asa 


linear combination of atomic orbitals age: 
ee ce, (V-19) 
v 


and Lf vs can be written as a linear combination of 
another basis set Xo then there are also linear combination 


GES the orbital St 
ee SC (V-20) 


where oe is a non-singular matrix. 


There are_three types of transformation: 
(a) Transformations which only mix atomic orbitals on 
the same atom with same n and 1. An important transforma- 
tion is the coordinate transformation. 
(b) Transformations which only mix atomic orbitals on 
the same atom with different 1. The resulting orbitals 
are known as atomic hybrid orbitals. 
Lc) Transformations which mix atomic orbitals on 
different atoms leading to a non-atomic basis set. 

For molecules of low symmetry coordinate invariance 
is an essential feature. However, invariance to hybridi- 


zation is less important. 
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Using the ZDO approximation in one frame, in 
another frame using the condition of equation (V-20) 


invariance can be proved as follows: 


(RPS) Leia ato a aren ear ETO oa (V-21) 
nvig HO VB AY o6 

If the zero differential condition is applicable to the 

old basis set, the integral (uv|Ac) will be 6 6 


uv Ao’ 
Therefore, 


(aB;ys) = } Sra Wet ky es (uu]AA) (V-22) 
u 


If we allow pure coordinate rotation only, then S does 
NOuycnange and “only Poy Po! Pp. rotate inside the sub- 


group. Therefore, 


S S 

= 1M 
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Of course, this was evident from the beginning. For the 


next step 


A = P, on B- (3 components) 


A’ B 
then 
(S,S,/y6) = } thy ty¢ (S,S, | AA) 
If we put 
Ts(ayp(p) ~ (Sa8qlPxpPxp) = (S,SalPypP yp) 


e (S,5,|PpP op) 


Then in the above equation 
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(S,S,1 76) = Ps (a) P(B) ) ie ie 
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SQmthate invariances is proved, 


ie 





7 a 
, iM 
= a 
megs cH, .,2= smwees ow 7) 
: A. Te" 


at A 


~&=<V eee Po | = 
_ / 4 fe ky (H) 3 A) S 
6 


| 
7 ' on : -. Ms a 
“I ieerenes, ‘Bk sonsizeval 
= oa 


Le H=Pay ASP, and if we put 


"p(a)e(B) = Px (ayPx (ay | Pac (B) Px (By) 





se (a) Poca) {Py (B) Py (By 


Px) Pe (ay | Pz (By? 2 (By 
ae aes bab Ma's el Sat, 3 xXeesoinanirola 


then (af|yé) = "p (A) P(B) 


rr 
~r 
ct 
15 
oe) 
ctr 
2 
RD 
cr 
~ 
~~ 
ct 
~ 
(or) 


aA (V-24) 


See 
P(A)p(B) %8 Y6 


Invariance is proved again. 

In the two center terms (u|V,|u) and (u[VR1v), 
monatomic differential overlap Ngads (utv) on atom A 
suggests that (ulv,|v) is equal to zero. According to 
the anvariance condition for all ay in A, the diagonal 

elements (u[V,,] 4) are the same. For the same reason 


as in’ [ we can write 
. AB 


(ulvV,1H) = Va, (v-25) 
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In the CNDO method the resonance integrals, Buve 
are evaluated semi~empirically. Here also the invariance 
condition has been satisfied by assuming that gp is 

vi 


proportional to the overlap integral. 


= 1 ° ° 
the eats (Bi 2H We ets (V~26) 


Thy Thy 


As overlap increases, the bonding capacity of the 
Overlap also increases. To satisfy the invariance 
condition it is required that for all atomic orbitals 
thespropomtionaliity factor between ao and ey is the 
Same. The constants are denoted by om 1 BY, and depend 
One the nature or atom itliands \: Benaket® may be prepared 
for each atom. In the present work we parametrize the 
theory in such a way that it encompasses these 


possibilities. 


Computational Method and Parametrization: 


We have developed an all-valence-electron computer 
program at the CNDO level of approximation. Our program 
Can wWiandle the snol lowing atomic orbitals andeaccmuc 


parameters: 
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Atomic Orbitals’: 
H 1s, 2p,, Dye eP 
Diy, ee, Dyn Oe 2S, 2Pyr eur 2p, 


Ne&peMg, Al, SiyePy Sy Cl 3979 3p_, oSel 3p, 


Atomic Parameters: 


Z (effective charge) 


ean eee ee 
S ee Py ee 
8B 
shy ore Tae te. 
FOr Mmoothnaan. ss, fLormila RSGeay) 
a a a a 2 
Sag pe Clepey erp ; 
x y Zz For Klondike formula a=(I-A) 
Cor ae i a. (Orbital Exponent) 
K (for Klondike formula) 
has 22 


The above 20 parameters can be assigned to each 
atom but kK, i, w are used as global parameters in the 
present work. 

Mostly we have used 1s orbitals for H and 2s, 2p. 
Nee and 2p, orbitals for carbon. We have also tested 
the effect of inclusion of higher orbitals such as 2p yr 
2BoF and 2P, orbitals on hydrogen in the case of benzene 
and we have compared the coefficients with those reported 


by Rerniaeen 
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Although agreement between these two sets is not 
So good, still we hope that inclusion of higher orbitals 
will give good results in the prediction of electronic 
transitiom leading to an excited state with considerable 
Rydberg character as Salahub and Sandorfy (33) did in the 
case of saturated paraffin. 

It should be obvious that the parametrization of 
the present computational program is a very flexible 
one. For a general molecular geometry the constraints 


placed upon the parameters due to invariance requirements 


may be classified as follows: 


(1) Both space (coordinate) invariance and hybridization 
invariance 
iy Se 
S he Py es 
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(2) Space invariance only 
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(3) Relaxation of hybridization invariance for hydrogen 
only 

For hydrogen case (2) 

For all others case (1) 
There are, however, a vast number of chemically important 
molecules which are of planar form with the o-7 distinction. 
These have been the SUD Ject= of applications of the 
immensely popular Pariser-Parr-Pople method. There is 
nor reason E£or=us “suddenly to abandon advantages due to 
the planarity in the all~valence-electron approximation. 
When the o-t distinction is possible, the constraints 
upon the parameters are 


(4) Both space and hybridization invariance for o part 
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(5) Space invariance only 
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For hydrogen case (5) 
For all others case (4) 


In the present work these possible flexibilities have been 


explored only to a limited extent. 
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Table (3) contains the values of CNDOR parameters 
along with CNDO/2 parameters for comparison. Aside from 
thesaddi tional global parameters, A and w, and the use 
Ofvuifferent @'setor 2s and 2p, there are no drastic 
changes in our parametrization from that of CNDO/2. 

| Table (4) contains the values of CNDOK parameters. 
A preliminary attempt has been made to parametrize CNDOK 
using k=0.4 which gives the Ohno's coulomb formula. 

The results of both methods have been presented 


in the next Chapter. 


719% 











Ve on 

i ™ ‘ 

2istemnieg AOI9D to astiisv ae andataog (£) sidst 
moxt obieA .nceliaygmos x02 Syevamarsq £\O0HD ditty ote 
; ue aan 

een of3 Snes \» bas |. ,exedomereq fadolp Lenoisiabs edd 


a 


otteeth on sxs ove .aS Bese ef 102 a'§ Seosezhie ae 
ss 
2 


7 
= 


, 
Syour> to 2st mon? noltesiazamstag 160 af eeGneaa 


v 5 


,aaelomeisg 4AKiMD to. sasulav oad snistno> HY sidst ; 

HO0MD exliayomsisg oF shem. need ead Janette yrsntml lexg A 

.sleero? dueafuoo e'omid eda asvinp doidw b ben geiege 
bodnevery esd sve aboddesa dgod to etiuaex-sat 


istgsd> txen edd FE 


Table (3) 


(Parameters in CNDO/2 and CNDOR) 


CNDO/2 CNDOR 

Bonding Parameter in ev (8) 

B,, (1s) chad eae 

B. (2s) 2180 20 
Bo (2P.) =B, (2p,,) =B, (2p) Zee Ly Raye) 
Core Parameter in ev (U) 

U,, (1s) DE aes) 17.06 

U, (2s) Comes 70.8399 
Dee Uakep) Uai2p_ ) Glee G25 203 


Average Ionization Potential and Electron Affinityaineey (a) 


a, (1s) e116 TLT6 
aq (2s) 14.051 14.051 
ao (2p) 5.572 5.572 


Orbital Exponent (Z) 


Sq (28) =C,, (2p) 1.625 1.8 


Core Charge (2Z) 
1 
Z,(1s) ub 


2 (28) =2), (2p) 4 4 


Global Parameter 
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Table (4) 


(Parameters in CNDOK) 


(a) Bonding parameter (8) in ev 


Bay 9.0 
Bo (25) =B, (2P,.)=B, (2p) =8, (2p,) fm 
(b) Core parameter (U) in ev 

U,, (1s) 13.00 

U, (2s) Died) 4 
Uet2p ) Walter) Uap, ) Ae is 
(c) One Center two electron integrals in ev 

(ir ee Oe Hees) Zeb 50 

Ciao rOor Cares) iL OPA RES: 

(e-Ay for C (zp) POR hs 
(d) Orbital exponent (7) 

Cy (1s) ‘ee 

Ga (2s) = 6, (2p) 18 


(e) Global parameter 


ho z= Au = -0.343 
Wo = Wr = 0.45a.u. 2 
KASS = 0.4 
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SCF iteration has been carried out in these 
computations until the difference in electronic energy 
between two successive iterations is less than Moe ye 
Convergence to this degree may be achieved within a few 
cycles. For improper bond lengths of molecules i.e. too 
ease or too small, some difficulty in convergence has 
been observed. Matrix diagonalization has been performed 
using the standard Jacobi method. ; 

Overlap integrals have been calculated according 
to the formulas of Mulliken et a1 (36) using the exponent 
values for carbon and Hydrogen noted in Tables (3) and (4). 

Two electron integrals have been calculated 
SBCCOLG INC LO ROeHean ts formula ‘13) in CNDOR and 
according to the Klondike formula in CNDOK. 

Bond angles .needed for the calculation of 


9) 


geometry have been taken from Wiberg's ‘4 paper. 
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Chapter yI 


RESULTS AND DISCUSSION 


A. CNDOR Method: 


Table .(5) .contains the calculated total energies 
for all the molecules studied along with the corresponding 
CNDO/2 results. It also contains calculated and observed 
heats of atomization. Table (6) contains the Oro caL 
energies calculated by CNDOR along with CNDO/2. Table (7) 
contains the variation of energy with poererens bond 
bengths # @Figures Gey ephdeent the variation of 
energy with bond lengths. The bond lengths with minimum 
energy from the figures are included in Table (8) and the 
corresponding energy is included in Table (5). Table (8) 
also contains the bond lengths of the molecules along 


with the experimental values for comparison. 
Heats of Atomization: 


From Table (5) it is clear that the energy 
calculated by the CNDO/2 method is much deeper than that 


found in the present work. 
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For comparison with experimental values we have 


Calculated heats of atomization in the following way: 
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= ' i] i] 1 1 
Beton Number of "C' atom x E(C) + Number of 'H' atom 
<i) 
Seon reeh stone ote! energy 


The heats of atomization calculated by the CNDOR 
method are in good agreement with the experimental values. 
On the other hand the agreement between observed heats of 
atomization and those calculated by CNDO/2 is exceedingly 
poor. The deviation of our results from the experimental 
is within 6% in all cases whereas the CNDO/2 deviation 


Varies, from LOS2. to 187s. 
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Orbital Energies: 
Table (6) contains the occupied orbital energies 
of the molecules studied. We have compared our results 


(35) 


with those of Kern rorebpenvzene. Lt 2s found that 


the orbital energies calculated by this method are in 


good agreement with those of Henn 2 


who performed an 
ab initio Hartree-Fock calculation of good accuracy. 

For comparison we have also included in the table the 
orbital energies calculated by the CNDO/2 method. Again 
our results are more reasonable, In particular, the 
ordering of the orbital energies obtained in the CNDO/2 
method does not agree with that of Kern whereas our 
results are in good agreement. Similar results Aes also 
found in the case of ethylene i.e. the orbital energies 
and the ordering calculated by the CNDOR are superior 

to those from CNDO/2 when compared to the results of 


(oe who performed SCF calculations 


Moskowitz and Harrison 
with Gaussian basis sets of different sizes. They have 
shown that the ordering does not change with different 


basis sets. The orbital energies of their largest 


basis set for CH, are included in Table (6).. 
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F. Occupied Orbital Energies for CoH (D3, symmetry) 
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Bond Lengths: 


On the otner hand, the bond lengths, Table (8) 
obtained in this work are in poorer agreement with 
experiment than those of the CNDO/2 method. Our 
values are still, in general, more reasonable than 
ae of Winerowees and we do not obtain the very 


deep total energies and the poor heats of atomization 


which result from the CNDO/2 method, 
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Variation of Energy at Different Bond Lengths 


(Bond length in A° and energy in a.u.) 
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B. CNDOK Method: 


Table (9) contains the energy values for all the 
molecules studied, Table (10) contains the energies 
calculated using different bond lengths. Figures are 
drawn to show the variation of energy with bond lengths 
(Figures 15-22). Bond lengths for which the energy is 
minimum, are included in Table (11) along with observed 
bond lengths, and the corresponding energy is included 
in Table (9) which also contains the calculated and 


observed heats of atomization. 
Heats of Atomization: 


To compare the heats of atomization with 
experimental values, they have been calculated in the 


following way: 
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AE tom Number ore Cl eatome xb (C) + 
Numbermor = Hivatoms: x 6 (4) 
AH tom = AE Jtom - Total Energy 


Although the heats of atomization calculated by 
the CNDOK method do not agree so well with experiment 
as in the previous method (CNDOR), Stia. tne, CNDOK 
method is in much better agreement than the CNDO/2 


method. 
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Bond Lengths: 


From Table (11) it is evident that bond length 
obtained by CNDOK method is not So good as in the case 
of CNDOR. But we should not be discouraged because 
this method is still in the stage of preliminary 
Sematietrization. We hope that after proper parametri- 


zation the prediction of bond lengths will be improved. 
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Variation of Energy with Different Bond Lengths in CNDOK 


(Bond length in A°, and Energy in a.u.) 
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Table (i) 


Comparison of Calculated and Observed Bond Lengths in CNDOK 


(a) 


Molecule Bond Bond Length Observed 
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Cn Conclusion. 


In view of the results presented here, we may 
conclude that good theoretical values for heats of 
formation can be obtained through the proper parametri- 
-Zation of a CNDO-level semi-empirical theory. The 
essential point we have brought forward is the importance 
of the balance between the one-electron energy (negative) 
and the two-electron energy (positive). Under the ZDO 
approximation, many terms in the energy formula are 
eliminated. If there arises any appreciable imbalance 
due to the effect of these neglected terms between the 
One- and two-electron parts, then we are bound to have 
e poor total—energy value.’ It is our opinion that this 
has been the case with most of the methods thus far 
proposed by various authors. We have attempted to 
maintain the balance by introducing two parameters, 

A and w, in the one-electron part of the diagonal 
energy matrix elements. 

As ea nea earlier, the use of the ZDO 
approximation was, to some extent, justified by 
introducing the notion of the orthogonalized atomic 
orbitals (for example Léwdin's orbitals) in the P-P-P 
theory but it seems to be difficult to maintain the 


same logic and subsequent justification in the all- 
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valence-electron theory. However we keep it as a 
necessary evil in order to make the calculations 
feasibilegtorslargesmolecu les. 

Although the bond lengths calculated by the 
_CNDOK method are not so good as those obtained by 
CNDO/2 method, the heats of atomization and» total— 
energies are much better than those of the CNDO/2 
method. Furthermore, we must emphasize that the 
CNDOK method does not give bond lengths which are so 
far afield as to render the method hopeless. On the 
Seed we believe that the CNDOK method represents 
a significant step toward solving Klopman's dilemma. 
Although in this method we can employ either of the 
three proposed ways (Roothaan's, Ohno's, or Mataga- 
Nishimoto's formula) to obtain coulomb integrals, we are 
not restricted to the use of the particular corresponding 
values of Si and Ks Ther flexibility, representedeby 
these parameters has yet to be exploited. Even more 
flexibility can be seen when we consider the effect of 
also varying and w along with the k's. The present 
work has shown that such a parametrization leads to 
a useful semi-empirical method. Considering that there 
still remains a considerable amount of the flexibility 
in the method we can be optimistic about the usefulness 
of the method in the simultaneous prediction of a Valery 


of molecular properties. 
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